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1. Introduction 

Very recently, experimental results were published [Tj on the measurement of forces on 
closed, asymmetric electromagnetic resonant cavities. These results add to previous 
claims in the same line by an independent researcher who did the first experiments [2]. 
Such claims were criticized by the scientific community mainly due to the proposed 
theoretical explanation, as Maxwell equations and Special Relativity clearly indicate 
that no force is possible without the emission of radiation from the cavity. On the 
other hand, it appears that General Relativity might allow for such kind of reaction- 
less propulsion, as exemplified and noted for the first time inp], where the low velocity 
limit of some warp drive spacetimes was analyzed. As indicated there, negative energy 
densities are required to accomplish that and, notably, some scalar fields present this 
possibility [4J. Of course, in order to have measurable effects similar to those reported, 
the coupling of the scalar field to matter or other fields acting as its source should be 
sufficiently strong, and this is precisely what has been proposed in [5] for the case of the 
electromagnetic field as source of the scalar field to explain discordant measurements of 
Newton gravitational constant. It is then only natural to wonder whether that theory 
(or a similar one) may account for the forces reported in resonant cavities. Of course, 
all this is highly speculative, and more prosaic explanations for these forces should be 
considered first. We proceed on the assumption that all spurious effects were accounted 
for, and on the belief that the possibility presented here is worth exploring. 

The theory put forward by Mbelek and Lachieze-Rey in [5] (see also [6]) represents 
a reduction to four dimensions of a Kaluza-Klein theory coupled to an external scalar 
V^, which in turn couples to matter. It is the source term of ip which allows for a 
possible strong coupling of the Kaluza-Klein scalar to other fields, in particular 
to the electromagnetic field. The theory was applied in cosmological[7j and galactic 
situations p] and, as mentioned, it was also used to investigate the possibility of the 
Earth's magnetic field influencing the measurements of Newton gravitational constant. 
In all these applications of the theory only its weak-field limit was used, and as this 
limit is similar for a wide range of theories, we employ in this work a rather general 
scalar-tensor theory, which incorporates the additional external scalar ip. 

In the next sections the equations of the mentioned scalar-tensor theory are derived 
from its proposed action, along with the equation of motion of neutral matter. Some 
axisymmetric electromagnetic modes of a truncated conical cavity are then presented 
and used as source in the weak-field approximation of the equations, previously obtained, 
to determine the force on the cavity. It is found that a coupling of the same magnitude 
as used in between the scalar and the electromagnetic field results in a correct 
magnitude and sign for the forces reported in asymmetric resonant cavities. As expected, 
the solution for the cavity presents negative energy densities (more precisely, it violates 
the weak energy condition [9]). The theory, however, does not seem to be completely 
satisfactory because in its linearized version it also predicts strong gravitational effects 
by the Earth's magnetic field, which are clearly not observed. A possible resolution of 
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this problem is considered in the last section. 



2. Scalar-tensor theory 

We will consider a scalar-tensor theory of the Brans-Dicke type [10] with inclusion of a 
Bekenstein's direct interaction of scalar and Maxwell fields[Il], and with an additional 
external scalar field ip minimally coupled to gravity, and universally coupled to matter, 
with action given by (SI units are used) 

levrGo J levrGo J (j) 
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+ \ j '^rnat [exp Q^^] dVt. (l) 

In order to have a non-dimensional scalar field (j) of values around unity, in expression 
(II]) the constant Go representing Newton gravitational constant is included, c is the 
velocity of light in vacuum, and Sq is the vacuum permittivity. Cmat is the lagrangian 
density of matter, which is assumed to couple to the scalar ip. The other symbols are 
also conventional, R is the Ricci scalar, and g the determinant of the metric tensor 
g^y. The Brans-Dicke parameters u (0) is considered a function of 0, as it usually 
results so in the reduction to four dimensions of multidimensional theories[T2]. The 
function A (0) in the term of the action of the electromagnetic field is of the type 
appearing in Bekenstein's theory and other effective theories[7j, it does not intervene in 
the weak field approximation ultimately employed, but is included for completeness. The 
electromagnetic tensor is Ffj_y = V^Ai, — VuA^, given in terms of the electromagnetic 
quadri-vector A^, with sources given by the quadri-current j'^. U and J are, respectively, 
the potential and source of the field ip. The source J contains contributions from the 
matter, electromagnetic field and the scalar 0. In order to build upon a concrete model 
we follow the proposal in|5] (convenient dimensional factors differing from those in[5] 
are employed here) 

J = /3mat (^, 0) ^T-'^* + ^EM (^, 0) ^^^F.^F^'^ + (^, 0) (2) 

where T™*^* is the trace of the energy-momentum tensor of matter, (note that this tensor 
is defined with respect to g^^^, not exp (/3'0) gf_iu), 

rjimat _ 2 6Cmat 

~ ^ Sg'^'^ ' 
and T'^ is the trace of the tensor 

= V^V,0 - V^V,<Pg^, + ^ (^V^0V.0 - ^ V>V^0(7;..) . 
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dQ 
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Variation of ([T]) with respect to (7'^^ results in {T^J^ is the usual electromagnetic 
energy tensor) 



SvrGn 



[A (0) T^^'' + T™/] + T^t 



(3) 



Variation with respect to gives 



0i? + 2a;V^V^^ 



- - — V>V,0 - 



) I 



dJ 



which can be rewritten, using the contraction of ([3]) with g^'^ to replace R, as 
{2uj + 3) VV^cj) 



du AnGpeo d\ SnGp , 

- -V 0V.0 - ^^<P^F,.F>^ + 



-rmat 



+ 



dJ 



-V^0, 



(4) 



where it was used that T^*^ = T/l^^g^"" = 0. 

The non-homogeneous Maxwell equations are obtained by varying ([T]) with respect 
to A^, 

V, {A (0) Fn = f^of- (5) 

with /io the vacuum permeability. 

The variation with respect to results in 



9-0 9-0 



(6) 



Having included Gq, it is understood that takes values around its vacuum 
expectation value (VEV) 0o = 1. The scalar -0 is also dimensionless and of VEV 

Finally, we consider the motion of neutral test particles, coupled to the scalar ip as 
indicated in ([1]), which is then obtained requiring that 



5 j mc\J exp (/30') g^^dx^^dx" = 0, 



to give 



Du^ 13 



Ds 



(7) 



It is important to mention that, in order to derive the previous equations, we have 
followed the prescription of not varying the trace of the energy-momentum tensors nor 
F^yF^^ in the source term ([2]), but only its coefficients /9's, as done in[5]. There is 
no clear reason for doing so, but on the one hand, inconsistent equations result if the 
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mentioned variations are included. On the other hand, the exact source term may not 
depend exphcitly on the tensors considered, and only after the equations are derived 
and substitutions made might it be expressible in terms of the say tensors. 



3. Weak-field approximation 

In the weak field approximation, for values of g^^ around r/^^^ taken as those of flat 
Minkowski space with signature (1,-1,-1,-1), so that g^i, = t]^,^ + h^^, we have 



with 



where 



h = rp^h^s = —rp^h-ys- 

The system ([3])-([S]) can then be written, to lowest order in h^^ and in the perturbations 
around the VEV's of and ip, as 

- V^'d^sKu = ^-^Tjr' + 2 {d,u<l> - V''d,s<pr^,.) , (8) 
with the Lorentz gauge 



d^hl = 0, (9) 

^0, (10) 



{2uo + 3) = ^T-»* - 

d.F^'' = /ioJ^ (11) 



V-o, (12) 



where ujq = c<;(0o)- It was used in these equations that, in order to recover the 
usual physics when the scalar fields are not excited, one must have A (0o) = 1) 
U {ipo) = J ii^o, 4>o) = 0. Also, as according to|5] the contribution form the energy- 
momentum tensor of the electromagnetic field through the source J is much larger than 
all its other contributions, the latter were neglected in the above equations. 

For slow moving neutral masses, the equation corresponds to the action of a 
specific force (per unit mass) (Latin indices correspond to the spatial coordinates) 



h = -j^ {hoo + hkk + 2/3^) + c^. (13) 
Introducing the D 'Alembertian operator 
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applying it to the force equation (IT^ . and using Eq. (jS]), one easily obtains 



mat 



I — \_f U / rpmat I rpmat \ 



+ 00 - 



2 9Xi V C2 

From Eqs. f lTOj) and f|T2|) . and retaining only the most important component T^"*, this 
expression can be conveniently recast as 

f = 

OXi 

where the "gravitational potential" x satisfies 



^J^rpmat , ^TtCq / 2 _ ^ T^mat 



+ + ■ (14) 



<t>o,i>o 



The first term corresponds to Newton gravity with only the matter source, while the 
second term, if one takes /3 = 0, corresponds to the matter contribution through the 
scalar 0, which is strongly constrained by Solar System tests, requiring large values of 
Uo- An interesting conclusion (not to be explored further here) is that the inclusion of 
the external scalar ip could thus allow ~ 1 if 

2uo + 3 

Note that the condition (1151) does not invalidate the conclusions in [S] as only the term 
dependent on the matter velocity of the "force" in the right hand side of ([7]) is used 
therein to explain the dynamics of rotating spiral galaxies. 

Making explicit the equation of the scalar 0, Eq. (ITU]) , with the expression of the 
source J, Eq. ([2]), one has 



mat 



[B' - E'/c') , (16) 



{2uo + 3) c2 

where it was used that, in terms of the modulus of the electric and magnetic vector 
fields, E and B, respectively, one has 

and where the contribution from itself as its source was not considered because, even 
if it is present, in the weak- field approximation one has 

T'^ = -300, 

and so its effect amounts to a redefinition of the rest of the coefficients in the equations 
for and ip. 
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It is also expected that the correction from the source J to the matter term in fll6l) 
is small, that is, 

dPmat 



^0 



< 1, 



if not, it should be included in the condition f|T5|) . Actually, condition f lTSj) refers to 
the motion of massive bodies, while the most stringent bounds on uq come from the 
propagation of electromagnetic waves near the Sun [13], not affected by the coupling of 
to matter. These bounds can be accomodated, always with uq ^ 1, ii 



1. 



(17) 



With all this, the contributions other than the matter to the potential x can then 
be obtained from ( IT^ as (we write x = Xmat + x') 



92 



8(3 



EM 



2uJo + 3 d(f) 



/3 



EM 



dip 



{B'-E'/c^). (18) 



Inp] it is argued that in order to explain discordant measurements oi G = Gq/ (j) a,s 
due to the generated by the Earth's magnetic field according to ( |T6|) . one must have 



SttGo^^o 



EM 



-(5.4 ±0.6) X 10" 



00, -00 



42 



(19) 



{2uJo + 3) c2 

while the value of /3 dJ/dip]^^^^ does not enter the equation of and is thus left 
unspecified. In the following we will evaluate the force predicted by ( ITSll for the resonant 
electromagnetic field in a conical cavity, assuming that the coefficient in the brackets in 
( 1T8|) can be estimated from the value ([T9|) alone. 



4. Normal modes in a conical cavity 



As done in [14] we consider a conical cavity with side walls corresponding to a truncated 
cone, with spherical sections as end caps. The cone axis is taken as the z direction, the 
lateral wall corresponds to the spherical angle 9 = 9q (half angle of the cone), and the 
spherical caps to the radii r = ri 2, with r2 > ri. 

The resonant modes correspond to standing electromagnetic waves satisfying the 
vector wave equation (F stands for either the electric field E, or the magnetic induction 



Bl 



1 d^Y 



- V^F = 0. 



The modes with rotational symmetry and B transverse to the z direction (called the 
TM modes) that satisfy this equation are (spherical coordinates are employed, with unit 
vectors e^, and e<^) (see [13] and references therein for details) 

B = -GkR{r)Q' {6) cos {ujt)e^, (20) 
Rir) 



E/c = C 



-n {n + l)Q (6) 
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+ 



R(r) 



+ R' (r) 



Q' (9) ee \ sin (ut) 



(21) 



where C is a global constant. The functions R and Q are defined as 

Q{e) =P„(cose), 

R (r) = i?+ (r) cos a + R^ (r) sin a, 

R^ (r) = 

where Pn is the Legendre polynomial of order n, Jm the Bessel function of the first 
kind of order m, and a and k constants to be determined along with the order n. 
By construction, the magnetic field satisfies the boundary condition of zero normal 
component at the metallic walls, while in order to have zero tangential components of 
the electric field at the walls, the order n of the Legendre polynomial must satisfy 

Pn{cOs9o)=0, 

the wavenumber k the condition 











V + 








r2 


r 


ri 




ri 


r 



r2 



and a 



tana 



R+ (ra) /r2 + R+ (r^) 



R- (ra) /r2 + R'- (^2) ' 
The resonant mode angular frequency is thus determined as a; = kc. 

There exists a complementary set of modes with E transverse to the z direction 
(TE modes), but for concreteness we study only the lowest frequency TM modes. 

An important parameter is the quality factor of the cavity, Qcavi for each mode. It 
is conventionally defined as 



Qcav — 



(W) 



(22) 



where u is the angular frequency of the mode, (U) is the temporal average of its 
electromagnetic energy, and (W) is the average dissipated power in the wall cavities. 
As the average electric energy is equal to the average magnetic energy in the cavity, 
and the loss power can obtained from the value of the magnetic field at the boundary, 
an explicit, practical expression of Qcav can be obtained in terms of solely the magnetic 
field as [15] 

_ 2 / {B') dV 

where the integrals are extended to the volume and the internal surface of the cavity, 
respectively, and S is the penetration length in the metal wall, of resistivity rj, 
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2r] 



(23) 
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From fl2UI) one can thus write 

2 j[Rir)Q'i9)fdV 

S J[R(r)Q'{9)fdS- ^ ' 

If the cavity is fed with an average electromagnetic power P, in the permanent 
regime one has (W) = P, and so, from (1201) and (1221) . 



J (52) dV _ C^k^ 



[R{r)Q' {e)fdV 



Qcav-P 



(25) 



which allows to determine the global constant C, given the fed average power and the 
characteristics of the cavity for the considered mode. 



5. Force on the cavity 

In the permanent regime of the established resonant mode, sustained against decay by 
a continuous power input P, the electromagnetic field (!20l) -(|2T]) corresponds to 



- E^jc^ = Fb (r, 9) cos^ (cot) - Fe (r, 9) sin^ {ut) 
= \{Fb- Fe) + \{Fb + Fe) cos {2ut) 



where 



FB{r,9) 
FE{r,9) 



C^e [R (r) Q' 



(26) 
(27) 



+ 



R(r) 



n{n + l)Q (9) 



R(r) 



+ R! (r) 



Q" (0) 



(28) 



From fll6p . expression (1261) then leads to a constant plus an harmonic in time 
contribution to (p, which, together with fl26l) in f|T8|) . result in x' also having a constant 
plus an harmonic part. The latter has a zero contribution to the time average of the 
force, and so we consider only the constant part, x'o; whose equation is, from (fTSjl . 

V\', = x{Fb-Fe), (29) 

where, using ( !T9l) . 



1 



\ 2wo + 3 

27rG'o£oV^o d(3EM 
2ujq + 3 d(f) 



8/3 



EM 



/3 



EM 



dip 



</>o,'/'o 



1.2 X 10' 



Am 



(30) 



Note that the magnetic field in the right-hand side of f|T8|) is the total field, which 
includes the contribution from the Earth's magnetic field. The latter, although of much 
smaller magnitude than that of the cavity cannot be neglected due to its large spatial 
scale. However, one has for the time average (denoted by (...)) 



B 



Earth 



B. 



cavity 



) 



B 



Earth 



B 



cavity 



> + 2B 



Earth 



(B 



cavity 



) 



B%arth + Fb/2, 



10 



since (BcavUy) = 0. In this way, the contribution from the magnetic fields of the Earth 
and of the cavity to the potential x'o can be separated, and that of the cavity alone is 
correctly described by (l29l) . 

Eq. ( 129|) is solved taking into account that its right-hand side is zero outside the 
cavity, so that, using the axial symmetry, the solution of Poisson equation (l29l) is 

Fn (r', 6') - (r', 6') 



X 

n 



^2 I ^/2 



X K 



2rr' cos ( 
4rr' sin 6 sin 6' 



r'2 sin 9' dr'de' 



(31) 



2rr' cos 

where K is the complete elliptic integral of the first kind, and the integral is extended 
to the interior of the cavity. Note that as the volume integral of the left-hand side of 
(!29|) is equal to zero, Vxo decays rapidly outside the cavity. 

Assuming a cavity with thin walls (but much thicker than the penetration depth 5, 
in order to the boundary conditions used to be correct) of mass surface density a, the 
force on the cavity is finally evaluated as 

F = -a y Vx'odS, (32) 

where the integral is extended to the internal surface of the cavity. Due to the axial 
symmetry the force has only a z component 

dx'o 



a 



a 



dz 



■dS 



cos 6 



dr 



/ Vxo ■ e^dS, 



dS, 



dx'o sin 6 dx'o 



de 



which is explicitly written as 



2no 



00 



cos 6 



dx'o sin 9 dx'o 



dr 



+ rl 



00 



cos 9 



dx'o 



' d9 
sin 9 dx'o 



sm 9d9 



r2 



dr 



+ sin 6*0 



r2 



COS 9 



dx'o 
dr 



■ d9 

sin 9 dx'o 
r d9 



sm9d9 



ri 



rdr. 



(33) 



00 



There are no details in the literature as to the precise dimensions of the cavities 
used in the experiments, so that an example roughly similar to the overall dimension 
reported and with the proportions observed in the published photographs will be used. 
Assuming a wall of thickness 1 mm, and a copper mass density of 8.9 x 10^ kg/m^, we 
have a = 8.9 kg/m?. 

We further consider the copper cavity to have ri = 18 cm, r2 = 36 cm, and 
^0 = 22°. For this cavity, the lowest TM mode corresponds to the order n = 5.75632 
of the Legendre polynomial, with a resonant frequency u = 1.05 GHz. For a resistivity 
rj = 1.72 X 10"*^ Q m the quality factor for this mode is Qcav = 3.13 x 10^. The next two 
TM modes have the same order n = 5.75632, and resonant frequencies u = 2.05 GHz 
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and V = 2.76 GHz, with quality factors Qcav = 3.11 x 10^ and Qcav = 5.24 x 10^, 
respectively. 

For an average power P = 1 kW, the constant C is evaluated for each mode using 
(!25|) . and ( 127|) and (128|) used in (!3T!) to obtain by numerical integration the values of 
Xq (r, 6^) needed in the numerical evaluation of (1331) . 

Note that, from fl25l) . the force on the cavity is proportional to the fed power, and 
to the quality factor Qcav 

For the lowest TM mode (z/ = 1.05 GHz) the value obtained is Fz = 7.7 N, while 
for the next two TM modes, with u = 2.05 GHz and u = 2.76 GHz, we obtained 
Fz = —1.4N and Fz = —0.9N, respectively. The values reported in[ij are not easy 
to compare with as the power of the microwave source is distributed over a rather 
wide range of frequencies, so that the actual power into the resonant mode is not 
precisely defined. Using a spectrum analysis of the power source the authors evaluate, 
for instance, that when Fz = — 0.3A^ the actual power into the resonant mode is 
P = 0.12 kW, which would correspond to Fz = -2.5 N at P = IkW. The last 
two modes considered are closer to the reported value of the resonance, u = 2.45 GHz, 
and give theoretical results with the correct sign and similar magnitude. As according 
to the model the force is proportional to the thickness of the wall, depends also on the 
precise geometry of the cavity (neither of them reported in the literature), and as the 
value (l30l) is only an estimation, since the contribution from (3 df3EM /dip\^^ cannot 
be ascertained independently, the results seem consistent with the measured force being 
due to the studied effect. 

Note that the lowest mode {u = 1.05 GHz) leads to a force much larger in magnitude 
and of opposite direction to that of the next two modes. This and other dependencies 
of the predicted force, as the proportionality to the cavity wall thickness (within certain 
limits as Vxo decays rapidly outside the cavity), can be explored experimentally with 
relative ease to test the theory. 

Finally, it is worth noting that the weak energy condition (WEC) [9] is violated for 
the cavity, as is the case in other models of propellant-less drive [3j. In effect, from ([8]), 
the WEC is written for the cavity 

{d^u<l> - ri^'d^5<pr]^.u) UW > 0, (34) 

for any time-like four-vector U^. By taking = (1, 0, 0, 0) one has the particular WEC 

V^^ > 0, 

which is seen from ( JT6i) and (|26ll to be violated at different times and regions inside the 
cavity. 

6. Discussion 

It was shown that the weak field approximation of a rather general scalar-tensor theory of 
gravity, which includes an additional scalar with strong coupling to the electromagnetic 
field, as proposed in|5], could account for the forces reported on asymmetric resonant 
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cavities. Although highly speculative, it is interesting that this was done using the 
same coupling coefficient adjusted by[5] to explain discordant measurements of Newton 
gravitational constant. It is also of interest that the inclusion of the external scalar ip 
can help to reconcile the Solar System tests with values of the Brans-Dicke parameter 
uj close to unity (see relation (ITTll ). The weakest part of the theory seems to be that 
there is no clear way of preventing large gravitational effects due to the magnetic field 
of the Earth, as predicted by Eq. f|T8|) . A possible solution can be sought in non-linear 
effects, such as those due to the second terms in the left-hand sides of (jl]) and ([6]). In 
effect, their inclusion would modify (fT4|) to 



□X = (Dx) 



original ~^ 2 



duj 



d''(f)d^(f) - I3d''(f)d^^ 



For the stationary case of the Earth's magnetic field one would then have 

1 




rpmat 
-^00 



If the terms V0 ■ V0 and V0 ■ Vip were to dominate over V^<; 
the equations (jl]) and dHD would result in 



(35) 



and V'^ipi respectively. 



duj 



V0- V0 



rj-imat _|_ 



dJ 



^0, 



<f'0,''l'0 



mat 



dJ_ 

dip 



00,-00 



which clearly cancel the terms in (l35|l leading to large values of the force. Put more 
plainly, and V^V cire sources of the potential x' and so situations where they 

are small or even zero would reduce or even nullify the gravitational effect of the 
electromagnetic field. Note that in the case of a static magnetic field outside its sources 
one can write B = V^', with V^\E' = 0, so it is possible that equations like (jl]) and ([6]) 
for the static case 

du 



(2wo + 3) V> + 



have the solutions V0 oc Vip oc V\E', and = V'^ip = 0, which is in general not 
possible for the case of the cavity, for which the constant part of — E'^/c^ cannot be 
written as V\E' ■ V\E'. 

If this is what happens in the case of the Earth's magnetic field, it would seem to 
invalidate the derivations in[5], where the solution with 7^ was used. However, 
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it can be shown that if du / d(t)\^^^^ / {2ujq + 3) ~ — 1 (the minus sign is necessary to 
obtain real solutions), both solutions are numerically similar. 

Note that if this is possible for a static magnetic field, it could not be possible for 
a static electric field outside its sources, which also satisfies E = V\l/, with V^\I' = 0, 
because the difference of sign would not allow real solutions. It is so expected that static 
magnetic fields show no unusual gravitational effects, while static electric fields and non- 
stationary electromagnetic fields do. Along these lines note finally that d'^(j)/dt'^is also 
a source of the potential x\ which would contribute in transient situations. 



Appendix: Modified Maxwell equations 



The generalized Maxwell equations, Eq. (|5]), are 



Since 



V, [A 
V. [A 



1 d 



A(0) 



one has 



d 



Keeping up to first order in the excitations of the scalars one has 
dPf"" , , , _ , „„„ d 



with 



dx'^ 



dx 



- {K<P - h/2) 



The equations for (p and h are at the same order of approximation (with only 
electromagnetic sources) 



^ttGo^^o , d(3EM 



0O>'/'O 



(2wo + 3) c2 
Dh = r^'^^D V = 600, 
so that, if one defines 9 = Aq0 — h/2, its equation is 



{B' - Eye') 



□9 



iirGoeo (Aq - 3) O/Sem 



(2wo + 3) c2 



(52 _ p2|^2^^ 



(36) 



<f'o,i>o 



while the Maxwell equations can be written in usual vector form as (neglecting the order 



one correction Xi 







60) to the four-current j'^) 
P 



V . E = ^ - V9 ■ E, 

£^0 
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dB 

VxE=-— , V-B = 0, 
ot 

^ ^ . IdE 1 dQ^ „ 

^^^ = ^°J + ^^ + ^^^-^®^^- 
These equations, together with Eq. f l36|) . form a complete set of for the electromagnetic 
field. 

It is interesting that in [5] the function A (0) obtained from reduction to four 

dimensions of the particular Kaluza-Klein theory employed is A (0) = 0^, so that 

Aq = 3, which according to Eq. (l36ll results in 9 = (the value of 9 is assumed 

zero at spatial infinity). Up to first order the Maxwell equations are so the usual ones 

in this case. For different theories, however, Aq 7^ 3, and the Maxwell equations can be 

modified. Note that even the case with no direct coupling, A (0) = 1, leads to modified 

Maxwell equations, which comes directly form the modification of the metric (relative 

to Minkowski space) due to the scalar field. 
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